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Abstract 

Direct photon production from a longitudinally expanding quark-gluon plasma 
(QGP) at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) 
energies is studied with a real-time kinetic description that is consistently incorpo- 
rated with hydrodynamics. Within Bjorken's hydrodynamical model, energy non- 
conserving (anti)quark bremsstrahlung q{q) — * q{q)^j and quark-antiquark annihi- 
lation qq — * 7 are shown to be the dominant nonequilibrium effects during the 
transient lifetime of the QGP. For central collisions we find a significant excess of 
direct photons in the range of transverse momentum 1 — 2 < p? < 5 GeV/c as 
compared to equilibrium results. The photon rapidity distribution exhibits a cen- 
tral plateau. The transverse momentum distribution at midrapidity falls off with a 
power law p^, u with 2.5<^<3asa consequence of these energy nonconserving 
processes, providing a distinct experimental nonequilibrium signature. The power 
law exponent v increases with the initial temperature of the QGP and hence with 
the total multiplicity rapidity distribution dN^/dy. 

Key words: Direct photon production, Quark-gluon plasma, Relativistic heavy ion 
collisions, Finite-temperature field theory 
PACS: ll.10.Wx, 12.38.Mh, 13.85.Qk, 25.75.-q 



1 Introduction 



The first observation of direct photon production in ultrarelativistic heavy 
ion collisions has been reported recently by the CERN WA98 collaboration in 
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208 Pb+ 208 Pb collisions at y/s = 158A GeV at the Super Proton Synchrotron 
(SPS) [1]. Most interestingly, a clear excess of direct photons above the back- 
ground photons predicted from hadronic decays is observed in the range of 
transverse momentum p T > 1.5 GeV/c in central collisions. As compared to 
proton-induced results at similar incident energy, the transverse momentum 
distribution of direct photons shows excess direct photon production in central 
collisions beyond that expected from proton-induced reactions. These findings 
indicate not only the experimental feasibility of using direct photons as a sig- 
nature of the long-sought quark-gluon plasma (QGP) [2] but also a deeper con- 
ceptual understanding of direct photon production in ultrarelativistic heavy 
ion collisions. 

Unlike many other new phases of matter created in the laboratory, the for- 
mation and evolution of the QGP in ultrarelativistic heavy ion collisions is 
inherently a nonequilibrium phenomenon [3,4]. Currently, it is theoretically 
accepted that parton-parton scatterings thermalize quarks and gluons on a 
time scale of about 1 fm/c after which the plasma undergoes hydrodynamic 
expansion and cools adiabatically down to the quark-hadron phase transition. 
If the transition is first order, quarks, gluons, and hadrons coexist in a mixed 
phase, which after hadronization evolves until freeze-out. Estimates based on 
energy deposited in the central collision region at the BNL Relativistic Heavy 
Ion Collider (RHIC) energies ^/s ~ 200A GeV suggest that the lifetime of the 
deconfined QGP phase is of order 10 fm/c with an overall freeze-out time of 
order 100 fm/c. Different types of signatures are proposed for each different 
phase. 

Of all the potential signatures of a QGP [5], direct photons and dileptons 
emitted by the QGP, i.e., electromagnetic probes, are free of hadronic final 
state interactions and can provide a clean signature of the early stages of a 
thermalized plasma of quarks and gluons. Therefore a substantial effort has 
been devoted to a theoretical assessment of the spectra of direct photons and 
dileptons emitted from the QGP [6-12]. 

The theoretical framework for studying direct photon production from a ther- 
malized QGP begins with an equilibrium calculation of the emission rate at 
finite temperature [8-12,14,15]. This rate is then combined with the hydrody- 
namic description to obtain the total yield of direct photons produced during 
the evolution of the QGP [8,9,11,16-18]. 

In a recent development, it was shown in Ref. [19] that direct photon pro- 
duction from a thermalized QGP with a finite lifetime, when calculated using 
a real-time kinetic approach, is strongly enhanced by energy nonconserving 
(anti)quark bremsstrahlung q(q) — > q(q)'~f and quark-antiquark annihilation 
qq — * 7. This nonequilibrium contribution has been missed by all the previous 
equilibrium rate calculations in the literature [8-10,12]. Most importantly, this 
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novel result poses a serious question on the applicability of the equilibrium 
rate calculations to direct photon production from a transient expanding QGP 
that is expected to be created in ultrarelativistic heavy ion collisions at RHIC 
and Large Hadron Collider (LHC) energies. 

The goals of this article, (i) We argue that the finite lifetime of a tran- 
sient QGP raises a conceptual inconsistency in the calculation of direct photon 
production via an equilibrium rate, (ii) We introduce a real-time kinetic de- 
scription which naturally accounts for the finite lifetime and nonequilibrium 
aspects of the QGP. (iii) This real-time kinetic approach is consistently com- 
bined with Bjorken's hydrodynamics [20,21] to obtain the direct photon yield 
for a QGP that is expected to be formed in central collisions at RHIC and 
LHC energies, (iv) We focus in particular on experimental signatures associ- 
ated with processes that would be forbidden by energy conservation in a QGP 
of infinite lifetime, (v) It is not our goal to assess photon production from the 
hadronic phase but to compare the real-time kinetic predictions for the QGP 
phase to those obtained from the equilibrium calculations. Furthermore we 
focus on extracting potential nonequilibrium signatures associated with the 
transient and nonequilibrium aspects of the QGP phase. 

Recently the equilibrium photon yield for a QGP with hydrodynamic expan- 
sion, along with the photon production from the hadronic phase has been used 
to estimate the total photon yield at CERN SPS energies [18]. This reference 
provides a thorough comparison between the yields from the QGP obtained 
from the equilibrium calculation, and that from the hadronic phase with the 
result that they are of comparable order. Instead our goal is to compare the 
nonequilibrium yield from the QGP to that obtained from the equilibrium for- 
mulation and we show a substantial enhancement from nonequilibrium effects 
indicating that the nonequilibrium yield from the QGP will stand out over 
that from the hadronic phase. 

Brief summary of the main results. We incorporate the real-time kinetic 
approach with the hydrodynamical evolution of the QGP to study direct pho- 
ton production for central collisions at RHIC (y/s ~ 200 A GeV) and LHC 
(y/s ~ 5500A GeV) energies. We find that during the finite QGP lifetime of 
order < 10 — 30 fm/c, energy nonconserving processes lead to a substantial 
enhancement in the total photon yield for 1 — 2 < pr < 5 GeV/c with a rapid- 
ity distribution that is fairly flat for \y\ < ?7 ccn , where |?7| < ?7 cen denotes the 
central rapidity region of the QGP. These processes also lead to a total photon 
yield at midrapidity that falls off with the transverse momentum as with 
2-5 < v < 3 in this region of transverse momentum. Numerical studies reveal 
that this exponent increases with the initial temperature of the thermalized 
QGP. This power law is a result of the transient and nonequilibrium aspects 
of the QGP and provides a distinct experimental nonequilibrium signature. 
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This article is organized as follows. In Sec. II we first review the usual S-matrix 
approach to direct photon production from a QGP. We then argue that this 
approach has shortcomings and is conceptually and physically incompatible 
with photon production from an expanding QGP with a finite lifetime. In 
Sec. Ill Bjorken's hydrodynamical model combined with the S-matrix rate 
calculation of photon production from the QGP is briefly summarized and the 
incompatibility of these two approaches is highlighted. In Sec. IV we introduce 
the real-time kinetic approach to photon production from a longitudinally ex- 
panding QGP which is consistently combined with Bjorken's hydrodynamics. 
This section contains our main results and we compare these to those obtained 
from the usual approach. Our conclusions are presented in Sec. V. In the Ap- 
pendix we show that the nonequilibrium photon production rate is related to 
the Fourier transform of the nonequilibrium expectation value of the quark 
electromagnetic current correlation function. 



2 S'-matrix approach and its shortcomings 

Production of direct photons from a QGP in thermal equilibrium has been 
studied extensively [6-10,12] because of its relevance as a clean probe of the 
QGP. Since photons interact electromagnetically with quarks their mean free 
path is much larger than the size of the QGP, hence they are not thermalized 
and leave the medium carrying direct information from the QGP. Because 
of the smallness of the electromagnetic coupling, the photon production rate 
is calculated to lowest order in the electromagnetic interaction. This rate is 
obtained from the S'-matrix calculation of the transition probability per unit 
time and is related to the Fourier transform of the thermal expectation value 
of the quark electromagnetic current correlation function [6,8,14,15], which in 
turn is determined by the imaginary part of the photon self-energy at finite 
temperature [8,14,15]. 

In order to highlight the shortcomings of this formulation and to establish 
contact with the real-time kinetic approach to photon production introduced 
in Sec. IV, we now review some important aspects of the calculation. 

Let us write the total Hamiltonian in the form 



where -£/qcd is the full QCD Hamiltonian, H 1 is the free photon Hamilto- 
nian, and Hi nt is the interaction Hamiltonian between quarks and photons 
with J M the quark electromagnetic current, the photon field, and e the 
electromagnetic coupling constant. 



H = Hq + Hi 




(2.1) 
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Consider that at some initial time ti the state \i) is an eigenstate of Hq with 
no photons. The transition amplitude at time tf to a final state \ f, 7a (p)) = 
\f) ® \l\(p)), again an eigenstate of H but with one photon of momentum p 
and polarization A, is up to an overall phase given by 

S(t f ,U) = (f, lx (p)\U(t f ,U)\t), (2.2) 

where U(tf, ti) is the time evolution operator in the interaction representation 



U{t f ,ti) =T exp -i H intJ (t)dt 

~ 1 - ie J*' dt J d 3 xJ?(x,t)A^ I ( y x,t) + £>(e 2 ), (2.3) 

where the subscript / stands for the interaction representation in terms of H . 
In the above expression we have approximated U(tf,ti) to first order in e, 
since we are interested in obtaining the probability of photon production to 
lowest order in the electromagnetic interaction. The usual S'-matrix element 
for the transition is obtained from the transition amplitude S(tf,ti) above in 
the limits ti — > — oo and tf — > oo 



Sfi = S(+oo, — oo) 

=-- 7 = y ^ y ^ dte ip ^e*(f\nx)\i), (2.4) 

where E = \p\ and P M = (£",p) are the energy and four-momentum of the 
photon, respectively, and is its polarization four-vector. Since the states 
\i) and |/) are eigenstates of the full QCD Hamiltonian Hqcb, the above S- 
matrix element is obtained to lowest order in the electromagnetic interaction, 
but to all orders in the strong interaction. We note that the S'-matrix element 
in effect is the amplitude for the transition between asymptotic states \i; in) — > 
I/, 7 A (p); out), i.e., S fi = (f, 7a (p); out|i; in), where |/,7 A (p);out) = |/;out) ® 
| 7a (p); out). Here, |7A(p*);out) is the asymptotic out state with one photon 
of polarization A and momentum p, and |i;in) (|/;out)) is the asymptotic in 
(out) state of the quarks and gluons. 

The rate of photon production per unit volume from a QGP in thermal equilib- 
rium at temperature T is obtained by squaring the ^-matrix element, summing 
over the final states, and averaging over the initial states with the thermal 
weight e~P Ei /Z((3), where (3 = 1/T, is the eigenvalue of H corresponding 
to the eigenstate \i), and Z((3) = J2i e~^ Ei is the partition function. Using the 
resolution of identity 1 = X)/ |/)(/|> the sum of final states leads to the elec- 
tromagnetic current correlation function. Upon using the translational invari- 
ance of this correlation function, the two space-time integrals lead to energy- 
momentum conservation multiplied by the space-time volume Q = V(tf — ti) 
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from the product of Dirac delta functions. The term tf — U — > oo is the usual 
interpretation of 2n5(0) in the square of the energy conserving delta functions. 

These steps lead to the following result for the photon production rate in 
thermal equilibrium [8,11] 



dN = 1 1 d 3 p v 0E 2 

where is the Fourier transform of the thermal expectation value of the 
current correlation function defined by 



= J d A xe ip - x (J»{$)J v {x))p. (2.6) 



In the expression above (•••)/? denotes the thermal expectation value. To 
lowest order in e 2 but to all orders in the strong interactions, W< v is related 
to the retarded photon self-energy 11^ by [14] 

e2 ^ = ^7F3Y Imn ?- ( 2 - 7 ) 
Thus, one obtains the (Lorentz boost) invariant photon production rate 

dN g^ R 1 



Kapusta et al. [9] and Baier et al. [10] calculated this photon production 
rate at one-loop order and showed that the processes that contribute to the 
energetic (E ^> T) photon production rate are gluon-to-photon Compton 
scattering off (anti)quark q(q)g — > q(q)l and quark-antiquark annihilation to 
photon and gluon qq — > g~f. Aurenche et al. [12] have recently extended the 
result to two-loop order. They found that the two-loop contributions to photon 
production rate arising from (anti)quark bremsstrahlung qq(g) — ■> qq{g)'~Y and 
quark-antiquark annihilation with scattering qqq(g) — > q(g)l are of the same 
order as those at one loop, and completely dominate the photon production 
rate at high photon energies [12]. For two light quark flavors (u and d), the 
invariant direct photon production rate to two-loop order from a QGP in 
thermal equilibrium at temperature T is given by [9,10,12] 



E 



dN 



d 3 p d A x 



5 aa« 



cq 



9 2tt 2 



T 2 e -E/T 



, 0.23E 16(J t -Jl) 
m — H 7, 



In 2 + — ) 

327 



(2.9) 
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where Jt ~ 1.11 and Jl ~ —1.07 [13]. We would like to emphasize that the 
thermal photon production rate (2.9) and in general the result (2.7) has two 
noteworthy features that are important to our discussions below: (i) The ther- 
mal rate is a static, time-independent quantity as a result of the equilibrium 
calculation, (ii) Emission of high energy photons is exponentially suppressed 
by the Boltzmann factor e~ E l T . 

We have reproduced the steps leading to Eq. (2.8), which is the expression for 
the photon production rate used in all calculations available in the literature, 
to highlight several important steps in its derivation in order to compare and 
contrast to the real-time analysis discussed below. 

(i) The initial states \i) are averaged with the thermal probability distribution 
at the initial time ti for quarks and gluons. In the usual calculation this 
initial time tj — > — oo and the initial state describes the photon vacuum and 
a thermal ensemble of quarks and gluons. 

(ii) The transition amplitude is obtained via the time evolution operator U(tf, ti) 
evolved up to a time tf and the transition amplitude is obtained by project- 
ing onto a state |/) at time tf, which in the calculation is taken tf — > oo. 
The sum over the final states leads to the electromagnetic current correla- 
tion function averaged over the initial states with the Boltzmann probability 
distribution, i.e., the thermal expectation value of the current correlation 
function. 

(iii) Taking tf — > oo and ti — > — oo and squaring the transition amplitude leads 
to energy conservation and an overall factor tf — ti. The rate (transition 
probability per unit time per unit volume) is finally obtained by dividing 
by V(tf—ti). The important point here is that taking the limit oitf—tj — > oo 
results in two important aspects: energy conservation and an overall factor of 
the time interval tf—ti. The resulting rate is independent of the time interval 
and only depends on the photon energy (and obviously the temperature). 

(iv) Main assumptions in the usual computations. In order to compare 
our methods and results with those obtained within the usual framework 
described above, it is important to highlight the main assumptions that are 
implicit in all previous calculations of photon production from a thermal- 
ized quark gluon plasma and that are explicitly displayed by the derivation 
above. Firstly, the initial state at ti (which in the usual calculation is taken 
to — oo) is taken to be a thermal equilibrium ensemble of quarks and gluons 
but the vacuum state for the physical transverse photons. Furthermore, the 
buildup of the photon population is neglected under the assumption that 
the mean free path of the photons is larger than the size of the plasma and 
the photons escape without rescattering. This assumption thus neglects the 
prompt photons produced during the pre-equilibrium stage. Indeed, Srivas- 
tava and Geiger [22] have studied direct photons from a pre-equilibrium 
stage via a parton cascade model that includes pQCD parton cross sec- 
tions and electromagnetic branching processes. The usual computation of 
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the prompt photon yield during the stage of a thermalized QGP assumes 
that these photons have left the system and the computation is therefore 
carried out to lowest order in a with an initial photon vacuum state. Obvi- 
ously keeping the pre-equilibrium photon population results in higher order 
corrections in a. In taking the final time tf to infinity in the S'-matrix el- 
ement the assumption is that the thermalized state is stationary, while in 
neglecting the buildup of the population the assumption is that the pho- 
tons leave the system without rescattering and the photon population never 
builds up. These assumptions lead to considering photon production only 
the lowest order in a, since the buildup of the photon population will nec- 
essarily imply higher order corrections. Although these main assumptions 
are seldom spelled out in detail, they underlie all previous calculations of 
the photon production from a thermalized quark gluon plasma, 
(v) The main reason that we delve on the specific steps of the usual computation 
and on the detailed analysis of the main assumptions is to emphasize that 
there is a conceptual limitation of this approach when applied to an expand- 
ing QGP of finite lifetime. The current theoretical understanding suggests 
that a thermalized QGP results from a pre-equilibrium partonic stage on a 
time scale of order 1 fm/c after the collision, hence for consistency one must 
choose ti ~ 1 fm/c. Furthermore, within the framework of hydrodynamic 
expansion, studied in detail below, the QGP expands and cools during a 
time scale of about 10 fm/c. Hence for consistency to study photons pro- 
duced by a quark gluon plasma in local thermal equilibrium one must set 
~ 10 fm/c. Hydrodynamic evolution is an initial value problem, indeed 
the state of the system is specified at an initial (proper) time surface (to 
be local thermodynamic equilibrium at a given initial temperature) and the 
hydrodynamic equations are evolved in time to either the hadronization or 
freeze-out surfaces. The calculation based on the S-matrix approach takes 
the time interval to infinity, extracts a time-independent rate and inputs 
this rate, assumed to be valid for every cell in the comoving fluid, in the 
hydrodynamic evolution during a finite lifetime. 

As stated in the introduction, however, the QGP produced in ultrarelativistic 
heavy ion collisions is intrinsically a transient and nonequilibrium state. It is 
therefore of phenomenological importance to study nonequilibrium effects on 
direct photon production from an expanding QGP with a finite lifetime with 
the goal of establishing potential experimental signatures. 

The current understanding of the QGP formation, equilibration, and subse- 
quent evolution through the quark-hadron (and chiral) phase transitions is 
summarized as follows. A pre-equilibrium stage dominated by parton-parton 
interactions and strong colored fields which gives rise to quark and gluon pro- 
duction on time scales < 1 fm/c [3]. The produced quarks and gluons thermal- 
ize via elastic collisions on time scales ~ 1 fm/c. Hydrodynamics is probably 
the most frequently used model to describe the evolution of the next stage 
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when quarks and gluons are in local thermal equilibrium (although perhaps 
not in chemical equilibrium) [20,21]. The hydrodynamical picture assumes lo- 
cal thermal equilibrium (LTE), a fluid form of the energy- momentum tensor 
and the existence of an equation of state for the QGP. The subsequent evo- 
lution of the QGP is uniquely determined by the hydrodynamical equations, 
which are formulated as an initial value problem with the initial conditions 
specified at the moment when the QGP reaches local thermal equilibrium, i.e., 
at an initial time ^ ~ 1 fm/c. The (adiabatic) expansion and cooling of the 
QGP is then followed to the transition temperature at which the equation of 
state is matched to that describing the mixed and hadronic phases [11,17,18]. 

Our main observation is that the usual computations based on S-matrix the- 
ory extract a time independent rate after taking the infinite time interval, 
which is then used in a calculation of the photon yield during a finite time 
hydrodynamic evolution. 

While we do not question the general validity of the results obtained via 
the S'-matrix approach, we here focus on the signature of processes available 
during the finite lifetime of the QGP and that would be forbidden by energy 
conservation in the infinite time limit. 



3 Bjorken's hydrodynamical model 

In order to highlight the conceptual limitation of the ^-matrix calculation for 
direct photon production for a transient QGP we now review the essential 
features of the hydrodynamic description [20,21]. For computational simplic- 
ity we work within Bjorken's hydrodynamical model of longitudinal expansion 
of the QGP [20], which is briefly summarized in this section. The main as- 
sumption in Bjorken's model is longitudinal Lorentz boost invariance in the 
central rapidity region of the QGP. This is motivated by the observation that 
the particle spectra for the secondaries produced in p+N and N+N collisions 
exhibit a central plateau in the rapidity space near midrapidity. For a longi- 
tudinally expanding QGP, it is convenient to introduce the proper time r and 
space-time rapidity i] variables defined by 



where t and z, respectively, are the time and spatial coordinate along the 
collision axis in the center of momentum (CM) frame. The transverse spatial 
coordinates will be denoted as xt, hence the space-time integration measure 
is given by 



\Jt 2 — z 2 , 7] 




1 , t + z 



(3.1) 



T = 



d A x = t air drj o9xt- 



(3.2) 
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Invariance under (local) longitudinal Lorentz boost implies that thermody- 
namic quantities are functions of r only and do not depend on rj. 

In Bjorken's scenario [20,21] the QGP reaches local thermal equilibrium at a 
temperature Tj at a proper time of order Tj ~ 1 fm/c after the maximum over- 
lap of the colliding nuclei. The initial conditions for hydrodynamical equations 
are therefore specified on a hypersurface of constant proper time Tj. The equa- 
tion of state for the locally thermalized QGP is taken to be that of the ultrarel- 
ativistic perfect radiation fluid (corresponding to massless quarks and gluons). 
The longitudinal expansion is described by the scaling ansatz v z = z/t, where 
v z is the collective fluid velocity of the hydrodynamical flow and describes 
free streaming of the fluid. Hence the space-time rapidity equals to the fluid 
rapidity. In terms of r and 77 the scaling ansatz implies that the four-velocity 
of a given fluid cell in the CM frame is given by = (cosh 7/, 0, 0, sinhr/) with 
u^u^ = 1. The conservation of total entropy leads to adiabatic expansion and 
cooling of the QGP according to the cooling law [20,21] 

T(r) = Tj i^j . (3.3) 

Hence, the QGP phase ends at a proper time ry = Tj(Tj/T c ) 3 , where T c ~ 160 
MeV is the quark-hadron transition temperature. At RHIC energies the initial 
thermalization temperature is estimated to be Tj ~ 200 — 300 MeV, which 
entails that the lifetime of the QGP phase is of order < 10 fm/c. 

Direct photon production in the S-matrix approach 

As discussed in detail above, the usual (^-matrix) calculation of direct pho- 
ton production from a hydrodynamically expanding QGP proceeds as fol- 
lows [9,11,16-18]. 

(i) First the rate of direct photon production is calculated within the S'-matrix 
framework described in the previous section, leading to Eq. (2.9) for the 
invariant rate. This expression for the rate describes the photon production 
rate in the local rest (LR) frame of a fluid cell in which the temperature is 
a function of the proper time of the fluid cell. The rate in the CM frame is 
obtained by a local Lorentz boost E — > P^u^ and the replacement T(t) — > 
T(r): 

dN 

d 2 p T dy t dr drj d 2 xx 

where pt and y are the transverse momentum and rapidity of the photon, 
respectively. 



E 



CM 



dN 
d 3 p d 4 x 



(3.4) 



LR. 
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(ii) The direct photon yield is now obtained by integrating the rate over the 
space-time history of the QGP, from the initial hypersurface of constant 
proper time to the final hypersurface of constant proper time t/ at which 
the phase transition occurs. This leads to the following form of the total 
direct photon yield in the CM frame for central collisions: 

[P%,T(r)], (3.5) 

LR 

where Ra is the radius of the nuclei and —i] ce n < i] < i] ccn denotes the central 
rapidity region in which Bjorken's hydrodynamical description is valid. The 
fact that the S'-matrix calculation for the rate results in a time-independent 
rate (a consequence of taking the infinite time interval, and hence assuming 
a stationary source as discussed above) determines that the only dependence 
of the rate in the LR frame on the proper time is through the temperature 
which is completely determined by the hydrodynamic expansion. 

It is at this stage that the conceptual incompatibility between the S-matrix 
calculation of the photon production rate and its use in the evaluation of the 
total photon yield from an expanding QGP of finite lifetime becomes manifest. 

The hydrodynamic evolution is treated as an initial value problem with a 
distribution of quarks and gluons in local thermal equilibrium on the initial 
hypersurface of constant proper time Tj ~ 1 fm/c. The subsequent evolution 
determines that the QGP is a transient state with a lifetime of order < 10 
fm/c. The direct photon yield is obtained by integrating the rate over this 
finite lifetime. The S-matrix calculation of the rate for a QGP in thermal 
equilibrium, on the other hand, implicitly assumes that 7$ — > — oo and 7/ — > oo 
as discussed above in detail. Therefore while the rate has been calculated by 
taking the time interval to infinity assuming a stationary source, it is integrated 
during a finite time interval to obtain the total yield. 

The question that we now address, which is the focus of this article, is the 
following: Is this conceptual incompatibility of physical relevance and if so what 
are the experimental observables? 

In order to answer this question and to assess the potential experimental 
signatures from nonequilibrium effects, we must depart from the S'-matrix 
formulation and provide a real-time calculation of the direct photon production 
rate based on nonequilibrium quantum field theory. 



dN 
d 2 p T dy 



CM 



f]ccn 



dN 



T I -vccn ^ E 'd 3 pd 4 x 
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4 Real-time kinetic approach 



Recently there has been substantial progress in the real-time approach to 
quantum kinetics in nonequilibrium quantum field theory [23-26]. The ad- 
vantage of this approach is that it allows to study the time evolution of the 
single (quasi)particle distribution functions for the relevant degrees of free- 
dom directly in real time as an initial value problem. The initial conditions 
for this initial value problem are specified in terms of the initial state of the 
plasma at the onset of the evolution. Furthermore, in the weak coupling limit 
a perturbative dynamical renormalization group method has been developed 
to resum directly in real time the perturbative expansions for the distribu- 
tion functions. The corresponding real-time dynamical renormalization group 
equations are the quantum kinetic (or quantum Boltzmann) equations [25,26]. 

This novel real-time kinetic approach has been applied to derive quantum 
kinetic equations in hot scalar [25] and Abelian gauge (scalar and spinor 
QED) [25,26] theories in a gauge invariant manner. 

Compared to the usual approach to quantum kinetics in which transition prob- 
abilities are computed using Fermi's golden rule and energy conservation, this 
real-time approach has the following noteworthy advantages: (i) It is capable 
of capturing energy nonconserving effects arising from the finite lifetime of the 
plasma, as completed collisions are not assumed a priori (ii) Because inverse 
time acts as an infrared cutoff, the real-time kinetic approach reveals clearly 
the infrared (threshold) singularities that lead to anomalous nonexponential 
relaxation, thus transcending the usual quasiparticle approximation [24-26]. 
(iii) Since both the real-time kinetic approach and hydrodynamics are formu- 
lated as an initial value problem, they can be incorporated consistently on 
the same footing. This last point proves very important in the hydrodynamic 
description of photon production. 

4-1 Nonexpanding QGP 

We begin our discussion with the calculation of the invariant photon produc- 
tion rate from a nonexpanding QGP. This calculation is relevant because the 
result is interpreted as the invariant photon production rate in the local rest 
frame of a fluid cell. The corresponding rate in the CM frame is obtained 
simply by a local Lorentz boost and the direct photon yield is obtained by 
integrating the rate over the space-time history of the QGP as explained in 
the previous section. 

Because of the Abelian nature of the electromagnetic interaction, we will work 
in a gauge invariant formulation in which physical observables (in the electro- 
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magnetic sector) are manifestly gauge invariant and the physical photon field 
is transverse (for details see Ref. [26]). 

The real-time kinetic approach begins with the time evolution of an initially 
prepared density matrix (see Refs. [24-26] and references therein). Consistent 
with the hydrodynamical initial value problem, we consider that at the initial 
time ti quarks and gluons are thermalized such that the initial state of the 
QGP is described by a thermal density matrix at a given initial temperature 

In order to compare our results with those obtained from the usual S'-matrix 
calculation, we will assume that photons are not present at the initial time, 
i.e., the initial state is the vacuum for physical transverse photons and neglect 
the build up of the photon population. We emphasize that this is not an ex- 
tra assumption in our treatment, but is one of the main assumptions in all 
previous calculations, we simply adopt this assumption for comparison. This 
assumption justifies a calculation of the yield to lowest order in a. While the 
usual approach assumes that the photons produced during the pre- equilibrium 
stage [22] had left the plasma without building up a population, this assump- 
tion can be relaxed by allowing an initial photon distribution and including 
the photon occupation number in the "gain" and "loss" terms in the kinetic 
equations [23-26] . Since the photons that were produced in the pre-equilibrium 
stage are a result of electromagnetic processes, the initial population will nec- 
essarily lead to higher order corrections in a. 

Thus consistently with all previous calculations and for comparison reasons 
we will assume no initial photon population and no population build up and 
compute the yield to lowest order in a. 

Therefore consistently with the S-matrix approach described above, the initial 
density matrix pi is taken to be of the form 

Pi = PQCD®|0 7 )<0 7 |, Pqcd = e-^ CD/T \ (4.1) 



We remark that the assumption that the initial density matrix is that of a 
thermalized system (after the strong interactions thermalize quarks and gluons 
on a time scale ~ 1 fm/c) underlies the program that studies the equilibrium 
properties of the QGP. This is our only assumption, i.e., that of a thermalized 
QGP at an initial time scale ti ~ 1 fm/c and is consistent with the general 
assumptions behind the equilibrium program. Of course this assumption is 
elevated to that of LTE, again consistent with the hydrodynamical description 
of an expanding QGP as an initial value problem. 
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At any time t later, the density matrix is given by 

p{t)=e~ iH{t - u) Pi e iH{t - u \ (4.2) 

with H the total (time independent) Hamiltonian given by Eq. (2.1). Since the 
initial density matrix is assumed to describe thermal equilibrium for quarks 
and gluons and therefore commutes with Hq defined in Eq. (2.1), it is straight- 
forward to find 

p{t) = e-^^Mijftr^Mje* (4.3) 

with U(t,ti) = e lH ot e ^H{t-ti) e -iH ti bging the time evolution operator in the 
interaction representation given by Eq. (2.3). 

We note that the upper and lower time limits in the unitary time evolution 
operator are a consequence of studying the time evolution of an initial state 
determined at ij up to a finite time t. This is the usual ingredient in the 
study of transition matrix elements in time-dependent quantum mechanics. 
In Refs. [25,26] we have previously provided a gauge invariant treatment of 
the real-time perturbative expansion for the Abelian sector which we use here 
to cast the real-time study of direct photon production directly in terms of 
gauge invariant observables (see also Ref. [19]). 

Taking 

PQCD = ^yE^ /Tl |j)0'l, (4.4) 

where \j) is the eigenstate of the full QCD Hamiltonian Hq C d with eigenvalue 
Ej and Z(Ti) = J2j e~ Ej / Ti is the initial partition function corresponding to 
Hqcd, we establish direct contact between the real-time kinetic approach and 
the S'-matrix approach discussed in the previous section. We haste to add, 
however, that unlike in the S-matrix formulation the initial and final times in 
the real-time kinetic approach will not be taken to =Foo, respectively. Keep- 
ing a finite interval of time as befits the description of an expanding QGP 
as a transient nonequilibrium state, as will be shown below, allows energy 
nonconserving processes that cannot be captured by the S'-matrix calculation. 

Since photons escape directly from the QGP without further interaction, it is 
adequate to treat them as asymptotic particles. The number operator N(p) 
that counts the number of photons of momentum p per unit phase space 
volume is defined by 

N(p) = j2a[(p)a x (p), (4.5) 

A=l 
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where a\(p) [a\{p)] is the annihilation (creation) operator that destroys (cre- 
ates) a photon of momentum p and transverse polarization A. These annihila- 
tion and creation operators can be written as usual in terms of the transverse 
component of the photon field and its conjugate momentum [19,26]. 

The number of photons per unit phase space volume at time t is given by 



1 1 d 3 pd 3 x 



= Tr 



Pit) Nip) 
PiNip,t) 



(4.6) 



where N(p, t) is the Heisenberg number operator. 

The invariant photon production rate is given by E dN{t)/d z pd A x and is ob- 
tained by using the Heisenberg equations of motion for the Heisenberg operator 
N(p,t) (for details, see Refs. [19,26]). A systematic framework to obtain the 
equation of motion is that of nonequilibrium quantum field theory in which 
the relevant correlation functions are obtained from a path integral along a 
contour in complex time [27]. A forward branch refers to the evolution for- 
ward in time and a backward branch refers to evolution backwards in time, 
and correspond, respectively, to the time evolution operator that pre- and 
post-multiplies the initial density matrix in Eq. (4.2). This method has been 
used previously in the treatment of quantum kinetic equations, and we refer 
the reader to Refs. [19,26] for details. The invariant photon production rate is 
given by [19,26] 



= hm (A - ie) V / A-Mi-k, 07 • A+(P, if) 

d 3 pd 4 x t>^t\dt> /^ 1 2(27r) 3 i (2vr) 3 v /V ' n TVF ' ; 

x^7(g,t))+c.c, (4.7) 

where k = p + q. In this equation, Nf is the number of quark flavors, e/ is the 
quark charge in units of the electromagnetic coupling constant e, At is the 
transverse component of the photon field, tpf is the (Abelian) gauge invariant 
quark fields, (• • • ) denotes the full nonequilibrium expectation value, and the 
superscripts "±" refer to fields defined in the forward (+) and backward (— ) 
time branches [27]. 

As mentioned above we assume that there are no photons initially and that 
those that are produced escape from the plasma without building up their 
population. Therefore the QGP is effectively treated as the vacuum for photons 
consistently with one of the main assumptions in all calculations performed 
in the literature. 
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Fig. 1. Feynman diagrams for the invariant photon production rate given by 
Eq. (4.7). Fig. (a) is in terms of the full self-energy of the photon to all orders 
in a s and to order a. Fig. (b) is the lowest order contribution, of order a. 

Consequently, the nonequilibrium expectation values on the right-hand side 
of Eq. (4.7) are computed perturbatively to order a and in principle to all 
orders in a s by using real-time Feynman rules and propagators. The photon 
propagators are the same as those in the vacuum. 

We shall further assume the weak coupling limit a <C a s < 1. Whereas the 
first limit is justified and is essential for the interpretation of electromagnetic 
signatures as clean probes of the QGP, the second limit can only be justified 
for very high temperatures, and its validity in the regime of interest can only 
be assumed so as to lead to a controlled perturbative expansion. 

In the appendix we show that invariant photon production rate given by 
Eq. (4.7) is related to the Fourier transform of the nonequilibrium expecta- 
tion value of the quark electromagnetic current correlation function to lowest 
order in a and to all orders in a s . The corresponding Feynman diagrams are 
shown in Fig. 1. 

We focus on the lowest order (one-loop) contribution which, as has been shown 
in Ref. [19], is missed by all the previous investigations which use the S- 
matrix approach thus assuming an infinite QGP lifetime. In using bare quark 
propagators we consider the quark momentum in the loop to be hard, i.e., 
q > Tj. Soft quark lines require hard thermal loop (HTL) resummed effective 
quark propagator [28] leading to higher order corrections. Indeed, the one- 
loop diagram with soft quark loop momentum is part of the higher order 
contribution of order aa s that has been calculated in Refs. [9,10]. 

To lowest order in perturbation theory, the result for two light quark flavors 
(u and d) reads [19] 
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with 



H(u) = / (0s{2[l - (P • fe)(p • q)]n(q)[l - n(k)] S(u + k - q) 

- (p • k)(p • q) n(q)n(k) 5{uj — k — q) >, (4.9) 



+ 

where q = \q\, k = \p + q\ and n(q) = \j\e q l Ti + 1] is the quark distribu- 
tion function at the initial time ij. The dependence of TZ(u) on the quark 
distribution at the initial time t{ is a consequence of the fact that in nonequi- 
librium quantum field theory the important ingredient is the initial density 
matrix [27,24-26], which consistent with hydrodynamics is taken to be ther- 
mal for quarks and gluons. Obviously TZ(cu) vanishes if there is no "medium" 
hence this contribution to photon production is a result of the presence of the 
medium, which exists as a thermal bath for t > U ~ 1 fm/c. 

At this stage we comment on how the assumption of neglecting the photon 
population can be relaxed. Allowing an initial as well as the buildup of the 
photon population leads to a full quantum kinetic equation of the form [24-26] 

-^n 7 (p, t) = r gain (p, t)[l + n 7 (p, t)] - T loss (p, t) n 7 (p, t), (4.10) 

where n 7 (p, t) = dN(t) / d 3 pd 3 x, r gain is given by the right hand side of Eq. (4.8) 
divided by E, and Ti oss is obtained from r gain through the replacement n <-> 
1 — n of the quark distribution functions (to the order considered) [24-26] . 
The population at the time of onset of local thermodynamic equilibrium is 
a result of electromagnetic processes during the pre-equilibrium stage [22], 
i.e., n 7 (p, t ) oc a, thus leading to higher order corrections in a. Therefore, 
to lowest order in a the neglect of the n 7 in the quantum kinetic equation is 
justified and leads to Eq. (4.8). 

In the case of a nonexpanding and thermalized QGP, the integral over t" in 
Eq. (4.8) can be carried out directly leading to 

E a ' r^^-^'-'Ji (4.1D 

d 3 pd 4 x (27r) 3 i-oo K ' tt(uj-E) ' V ; 



A detailed analysis [19,26] of TZ(u) shows that the first delta function 5(u + 
k — q) with support below the light cone {oo 2 < E 2 ) corresponds to the Lan- 
dau damping cut, and the second delta function S(u — k — q) with sup- 
port above the light cone corresponds to the usual two-particle cut. Fur- 
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thermore, TZ(u) has a clear physical interpretation in terms of the follow- 
ing energy nonconserving photon production processes: the first term de- 
scribes (anti)quark bremsstrahlung q(q) — * g(g)7, and the second term de- 
scribes quark-antiquark annihilation to photon qq —> 7. It has been shown in 
Ref. [19] that the dominant process for production of energetic photons is the 
(anti)quark bremsstrahlung. As explained in Ref. [19], in writing 1Z(uj), we 
have ignored the term corresponding to the "vacuum" process — > qcpf. The 
"vacuum" term, which persists for an infinitely long time, results in an energy 
conserving delta function that obviously vanishes. 

At this stage we can make contact with the S'-matrix calculation and highlight 
the importance of the finite-time, nonequilibrium analysis. If, as is implicit in 
the S-matrix calculation, the QGP is assumed in thermal equilibrium and 
with an infinite lifetime (entailing that ti — > —00 and tf — > 00), then we can 
take the infinite time limit ti — > —00 in the argument of the sine function in 
Eq. (4.11) and use the approximation 



This is the assumption of completed collisions that is invoked in time-dependent 
perturbation theory leading to Fermi's golden rule and energy conservation. 
The delta function 5(oo — E) is a manifestation of energy conservation for each 
completed collision. Under this assumption one finds a time-independent pho- 
ton production rate proportional to 71(E), provided that the latter is finite. 
In the present situation, however, the delta functions in TZ(uj) cannot be satis- 
fied on the photon mass shell. Therefore, under the assumption of completed 
collisions the lowest order energy nonconserving contribution to the photon 
production rate simply vanishes due to kinematics. Therefore this lowest order 
contribution is absent (by energy conservation) in the S'-matrix calculation, 
but is present at any finite time. 

The relevant question to ask is how this finite-time contribution of order a 
compares to the higher order S-matrix contribution to the photon yield. 

The photon yield (per unit volume) is obtained by integrating the rate over 
the lifetime of the QGP. Using Eq. (4.11), one obtains 



where t f - U ~ lOfm/c is the lifetime of the QGP. For T; ~ 200 MeV, the 
following important results have been shown in Ref. [19]. 

(i) Direct photon production given by Eq. (4.13) features a power law spectrum 
for E ^> Ti as a consequence of the energy nonconserving photon production 



sin[(cj — E)t\ t^. 
7r(u - E) 



oo 



5(uj-E). 



(4.12) 




(4.13) 
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processes: (anti)quark bremsstrahlung q(q) — > q(q)j and quark- ant iquark 
annihilation qq — > 7. 

(ii) The direct photon yield arising form this nonequilibrium finite-lifetime ef- 
fect in the energy range E > 2 GeV dominates over that obtained from 
equilibrium rate calculations given by Eq. (2.9) during the QGP lifetime. 

(iii) Furthermore the general analysis presented in Refs. [25,19] shows that the 
general expression for the invariant rate given by Eq. (4.13) in terms of 
spectral representations is correct to all orders. In particular to lowest or- 
der in a but to all orders in a s the spectral density lZ(u) = a[1Z^\uo) + 
a s hi(l/a s )TZ^(uj) + •••], all terms enter with the same time dependence. 
In particular 1Z^(uj) is given by Eq. (4.9) and TZ^(u = E) is given by 
Eq. (2.9). In the limit tf — U — > 00 the terms with support at uo = E result 
in a photon yield that grows linearly with time, while the terms for which 
7Z(E) = will grow in time much slower. In the infinite time limit those 
terms that do not vanish at u = E lead to photon yield that grows linearly 
with time and hence a finite time-independent rate, while those that vanish 
will be subleading. However, for a finite time interval which terms in the 
perturbative expansion dominate will depend on the competition between 
the order in perturbation theory and the length of the time interval. 

While these results in the nonexpanding case revealed the importance of the 
nonequilibrium and finite-lifetime aspects, the most experimentally relevant 
case to study is that of an expanding QGP. A constant rate obtained from the 
S'-matrix calculation leads to a yield that will grow in time linearly, while as 
analyzed in Ref. [19] the lowest order nonequilibrium yield grows much slower 
(logarithmically). 

Therefore the relevant question is which contribution is dominant during the 
lifetime of the QGP . This question can only be answered by a detailed nu- 
merical study of both yields which is performed below. 



4-2 Longitudinally expanding QGP 



As a prelude to photon production from an expanding QGP, and according 
with the usual approach [17,18], we first focus on the invariant photon pro- 
duction rate from each individual fluid cell of the QGP. Since the proper time 
equals the local time in the local rest frame of any fluid cell, we can follow the 
same real-time nonequilibrium analysis presented in the proceeding subsection 
to calculate the nonequilibrium invariant photon production rate in the local 
rest frame of each fluid cell [17,18]. We obtain [19,26] 



19 



E 



dN(r) 
d 3 p d 4 x 



LR 




xf4(?»/ (g,r)) + c.c. 



(4.14) 



As before, the nonequilibrium expectation values on the right-hand side of 
Eq. (4.14) is computed perturbatively to order a and in principle to all orders 
in a s by using real-time Feynman rules and propagators. To lowest order in 
perturbation theory, the result for two light quark flavors (u and d) reads 



where Tt(u) is the same as that given in Eq. (4.9), but now with n(q) = 
\j\e q l Ti + 1] being the quark distribution function at initial proper time Tj, at 
which the QGP reaches local thermal equilibrium. 

In principle in the case of an expanding QGP under consideration, the photon 
production rate given by Eq. (4.15) has to be supplemented by kinetic equa- 
tions that describe the evolution of the quark and gluon distribution functions 
so as to setup a closed set of coupled equations. However the assumption of 
the validity of (ideal) hydrodynamics entails that the quarks and gluons form 
a perfectly coupled fluid. This in turn implies that the mean free paths of 
the quarks and gluons are much shorter than the typical wavelengths and the 
relaxation time scales are much shorter than the typical time scales, i.e., the 
quark and gluon distribution functions adjust to local thermal equilibrium in- 
stantaneously. Thus the assumption of the validity of (ideal) hydrodynamics 
bypasses the necessity of the coupled kinetic equations: quarks and gluons 
are in local thermal equilibrium at all times. Therefore, within the frame- 
work of hydrodynamics we obtain the invariant photon production rate by 
directly replacing the initial quark distribution function n(q) in Eq. (4.9) by 
the "updated" distribution function n[q, T{t")\ = l/[e' 3//T(r "* ) + 1] at proper 
time t" > Ti, where T{r") is determined by the cooling law Eq. (3.3). The 
nonequilibrium invariant photon production rate that is consistent with the 
underlying hydrodynamics is then given by 



E 



dN(r) 
d 3 p d 4 x 



LR 




dN 



[r, E, T{r)\ = 7 ^ V J T dr" f °° ^ K[u>, T{r")\ 



E 



d 3 p d 4 x 



LR 



x cos[(uj - E)(t - t")}, 



(4.16) 



with 
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3 / (27r) 3 l"L"*" ~ (p ' k)(P' <i)\ n [Qi T{r)\n[k, T{t)\ 
x S(u + k - q) + [l + (p ■ k) (p ■ q)]n[q, T(r)]n[k, T(r)\ 
x<J(w-Jfe-g)j, (4.17) 

where n[g, T(r)] = l/[e 9//T(T ) + 1] and n = 1 — n. The momentum g integrals 
in Eq. (4.17) are calculated in the LR frame and hence are equivalent to those 
of the nonexpanding case above. 



Before proceeding further, we emphasize two noteworthy features of the non- 
equilibrium invariant photon production rate: (i) The photon production pro- 
cesses do not conserve energy, (ii) The rate depends on (proper) time not only 
implicitly through the local temperature but also explicitly. Furthermore, this 
explicit (proper) time dependence is non-Markovian as clearly displayed in 
Eq. (4.16). Whereas these two features seem to be rather uncommon within 
the framework of the usual ^-matrix approach to transport in heavy ion col- 
lisions [29], they are not unusual in quantum kinetics of nonrelativistic many- 
body systems [23,30]. Indeed, energy nonconserving transitions and memory 
effects which cannot be explained by usual (energy conserving) Boltzmann 
kinetics have been observed recently in the ultrafast spectroscopy of semicon- 
ductors that are optically excited by a femtosecond laser pulse [31]. 

The real-time kinetic approach when incorporated consistently with hydro- 
dynamics reveals clearly that photon production from an expanding QGP is 
inherently a nonequilibrium quantum effect associated with the finite lifetime 
of the QGP. 

At RHIC and LHC energies the quark distribution function n[q, T(r)] depends 
on the proper time r very weakly through the temperature T(r) within the 
lifetime of the QGP phase, hence a Markovian approximation (MA) in which 
the temperature in TZ[u,T(t")] is taken at the upper limit of the integral is 
reasonable and hence the memory kernel may be simplified. In this Markovian 
approximation, 1Z[cu,T(t")] in Eq. (4.16) is replaced by 1Z[cu,T(t)] and taken 
outside of the r"- integral. Thus Eq. (4.16) becomes 



E 



dN 



d 3 p d A x 



MA 



[t,E,T(t)] 



LR 



qIuK[lu,T(t)] 



(2tt) 3 
sin[(u -E)(r - Ti)] 



x 



tt(uj - E) 



(4.18) 



A computational advantage of this Markovian nonequilibrium production rate 
is that it provides the "updated" quark distribution functions locally in time. 
Physically, the motivation for this approximation is that the most important 
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aspect of the nonequilibrium effect is the nonconservation of energy (i.e., off- 
shellness) originated in the finite lifetime of the QGP, a feature that is missed 
by the usual S'-matrix calculation, while the proper time variation of the tem- 
perature is a secondary effect and accounted for in the S'-matrix approach. 

It is worth noting that a connection with the Boltzmann approximation can be 
obtained by assuming completed collisions, i.e., taking the limit Tj — > — oo in 
the argument of the sine function in Eq. (4.18) and using the approximation 
given by Eq. (4.12). Consequently, the lowest order nonequilibrium photon 
production rate vanishes in the Boltzmann approximation due to kinematics. 
This highlights, once again, that the usual approach to photon production 
outlined in Sec. 2 and used in the literature corresponds to the Boltzmann 
approximation and therefore fails to capture the energy nonconserving and 
memory effects that occur during the transient stage of evolution of the QGP. 

We are now in a position to calculate the direct photon yield from a lon- 
gitudinally expanding QGP. In the CM frame the invariant production rate 
for photons of four-momentum P M from a fluid cell with four-velocity can 
be obtained from Eq. (4.18) through the replacement E — > P^u^. In terms 
of the photon transverse momentum px and rapidity y, one finds P^u^ = 
p T cosh(y — rj). The photon yield is obtained by integrating the nonequilib- 
rium rate over the space-time evolution of the expanding QGP. Assuming a 
central collision of identical nuclei, in the Markovian approximation we find 
the invariant nonequilibrium photon yield to lowest order in perturbation the- 
ory to be given by 

MA 

[r,P%,T(r)], (4.19) 

LR 

where Ra is the radius of the nuclei and —i] ccn < rj < r] cen denotes the central 
rapidity region within which Bjorken's hydrodynamical model is valid. 

As remarked above, at RHIC and LHC energies the quark distribution function 
n[q, T(t)] depends on the proper time r very weakly through the temperature 
T(t) within the lifetime of the QGP phase, therefore the resultant photon yield 
is expected to qualitatively resemble the nonequilibrium photon yield from 
a nonexpanding QGP studied in Ref. [19]. This will be numerically verified 
below. We note that the expanding and nonexpanding cases differ mainly by 
the Jacobian r in the r integral in Eq. (4.19) that accounts for the longitudinal 
expansion of the QGP, and by the replacement E — > P^u^ in the argument of 
the invariant photon production rate that accounts for the shift of the photon 
energy under local Lorentz boost. 



dN 
d 2 pT dy 



= 7rR A 



CM 



VcCn d E ^ N 
Vccn d 3 p d 4 x 
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Au+Au @ RHIC 




p T (GeV/c) 

Fig. 2. Comparison of nonequilibrium (solid) and equilibrium (dotted) photon yields 
at midrapidity (y = 0) from a longitudinally expanding QGP at RHIC energies with 
initial conditions given by Tj = 1 fm/c and T; L = 200 (top), 300 (bottom) MeV. 

4-3 Numerical analysis: central collisions at RHIC and LHC energies 

We now perform a numerical analysis of the nonequilibrium photon yield and 
compare the results to the equilibrium one obtained from higher order equilib- 
rium rate calculations given by Eq. (2.9). The nonequilibrium photon yield in 
the Markovian approximation given by Eq. (4.19) contains a four-dimensional 
integral that is performed numerically for the values of parameters of relevance 
at RHIC and LHC energies. 

For central 197 Au+ 197 Au collisions at RHIC energies y/s ~ 200 A GeV, we 
take Ra — 1.2 A 1 ^ fm 7 fm [15] and r] cen = 2. The initial thermalization 
time is taken to be 7* = 1 fm/c [18,20], the final proper time t/ is determined 
when the critical temperature for the quark-hadron transition is reached at 
T(tj) ~ 160 MeV and is obtained from the cooling law given by Eq. (3.3) for 
a given initial temperature Tj at proper time Tj. 

The nonequilibrium photon yield at midrapidity (y = 0) in the range of trans- 
verse momentum Tj < p T < 5 GeV/c is shown on a log-log plot in Fig. 2 
for initial temperatures Tj = 200 and 300 MeV. For comparison we also plot 
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Fig. 3. Rapidity distribution of the nonequilibrium photon yield at px = 1 (solid), 2 
(dotted), and 3 (dashed) GeV for a longitudinal expanding QGP at RHIC energies 
with initial conditions given by Tj = lfm/c and Tj = 200 (top), 300 (bottom) MeV. 
The distribution is symmetric at y = 0. 

the corresponding equilibrium yield obtained by integrating Eq. (2.9) with the 
transformation to the CM frame as specified by Eq. (3.5), and using the value 
of a s [17,18,32] 

a * [T(T)] = (33-2^)M8T(r)/T c ] ' (420) 

with Nf = 2. Furthermore, Figure 3 shows the rapidity distribution of the 
nonequilibrium photon yield at different values of p T for Tj = 200 and 300 
MeV. Several noteworthy features are gleaned from these figures: 

(i) Whereas the equilibrium yield dominates the total yield for small px, the 
nonequilibrium yield becomes significantly dominant in the range pr > 
1.0 — 1.5 GeV/c. Perhaps coincidentally, this is the range in which the 
CERN WA98 data for central Pb+Pb collisions at SPS energies shows a 
distinct excess [1]. 

(ii) While the equilibrium yield leads to a transverse momentum distribution 
that falls off approximately with the Boltzmann factor e~ PT ^ Ti , the nonequi- 
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librium yield is not Boltzmann suppressed and falls off algebraically. It is 
found numerically that for ^> Tj and within the region 1 < p? < 5 
GeV/c the nonequilibrium yield at midrapidity falls off with a power law 
Pt" with v ~ 2.47 and 2.77 for T t = 200 and 300 MeV, respectively. This is 
a remarkable consequence of the fact that bremsstrahlung is the dominant 
process. As discussed in detail in Ref. [19] at large energies the dominant 
process is bremsstrahlung corresponding to the contribution from the term 
n(q)[l — n(k)] with k = \p + q\ in Eq. (4.9). For large photon energy p the 
important contribution, which is not exponentially suppressed arises from 
the small q region. 

(iii) A central plateau in the range of rapidity y < 2 is seen clearly. The rapidity 
distribution begins to bend down when y ~ i] cen , i.e., when the photon 
rapidity probes the fragmentation region. 

(iv) The numerical analysis of the expanding case reveals features very similar 
to those found in the nonexpanding case studied in Ref. [19], where we 
estimate that for px > 1.5 GeV/c the higher order equilibrium contribution 
to the direct photon yield (which grows linearly in time) becomes of the 
same order as the lowest order nonequilibrium contribution (which grows at 
most logarithmically in time) only if the lifetime of the QGP is of the order 
longer than 1000 fm/c. Therefore, we conclude that these nonequilibrium 
effects dominate during the lifetime of the QGP in a realistic ultrarelativistic 
heavy ion collision experiment. 

Similar numerical analysis can be performed for central 208 Pb+ 208 Pb colli- 
sions at higher LHC energies y/s ~ 5500A GeV, for which we take Ra ~ 7 fm, 
f] cen = 5 and r, = 1 fm/c [18,20]. The comparison of nonequilibrium and equi- 
librium photon yields at midrapidity is displayed in Fig. 4 for Tj = 450 and 
500 MeV. The dominance of the nonequilibrium yield at high remains but 
now at higher transverse momentum pt > 2 GeV. This can be understood as a 
consequence of the longer QGP lifetime resulting from higher initial tempera- 
ture at LHC energies. Furthermore, it is found numerically that for 1 < pt < 5 
GeV/c the nonequilibrium yield at midrapidity falls off with a power law p^ v 
with v ~ 2.52 and 2.56 for T t = 450 and 500 MeV, respectively. 

These results indicate a clear manifestation of the nonequilibrium aspects of 
direct photon production associated with a transient QGP of finite lifetime. 
The most experimentally accessible signal of the nonequilibrium effects re- 
vealed by this analysis is the power law falloff of the transverse momentum 
distribution for direct photons in the range 1 < pt < 5 GeV/c with an expo- 
nent 2.5 < v < 3 for temperatures expected at RHIC and LHC energies. Our 
numerical studies reveal that this exponent increases with initial temperature 
and therefore with the initial energy density of the QGP and the total multi- 
plicity rapidity distribution dN n /dy [20]. This could be a clean experimental 
nonequilibrium signature of a transient QGP since the photon distribution 
from the hadronic gas is expected to feature a Boltzmann type exponential 
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Fig. 4. Comparison of nonequilibrium (solid) and equilibrium (dotted) photon yields 
at midrapidity (y = 0) from a longitudinally expanding QGP at LHC energies with 
initial conditions given by Tj = 1 fm/c and T; L = 450 (top), 500 (bottom) MeV. 

suppression for p T 3> Tj [9,16-18]. 



5 Conclusions and discussions 



Our goal in this article is to search for clear experimental signatures of direct 
photons associated with nonequilibrium aspects of the transient quark-gluon 
plasma created in ultrarelativistic heavy ion collisions at RHIC and LHC en- 
ergies. 

We argue that the usual ^-matrix approach to direct photon production from 
an expanding nonequilibrium QGP has conceptual limitations. Instead, we 
introduce a real-time kinetic approach that allows a consistent treatment of 
photon production from a transient nonequilibrium state of finite lifetime. 

We focus on obtaining the direct photon yield from a thermalized QGP un- 
dergoing Bjorken's hydrodynamical expansion in central Au+Au collisions at 
RHIC (y/s ~ 200 A GeV) and LHC (y/s ~ 5500 A GeV) energies. The life- 
time of a QGP for these collisions is of order < 10 — 30 fm/c. We find that 
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(anti)quark bremsstrahlung q(q) — > 9(^)7 and quark-antiquark annihilation 
gg — > 7, 6ot/i 0/ lowest order in a, dominate during such short time scales, 
with bremsstrahlung dominating for ^> Tj. The contribution from these 
processes is a consequence of the transient nature and the finite lifetime of the 
QGP. As compared to the equilibrium rate calculations, these processes lead 
to a substantial enhancement in direct photon production for 1 — 2 < px J$ 
5 GeV/c near midrapidity. In striking contrast with the equilibrium calcu- 
lation that predicts an exponential suppression of the transverse momentum 
distribution for p T ^> Tj (Tj is the initial temperature), the nonequilibrium 
processes lead to a power law behavior instead. We find that at RHIC and LHC 
energies the direct photon transverse momentum distribution near midrapid- 
ity is of the form p^? with 2.5<z/<3forl<py<5 GeV/c and that photon 
rapidity distribution (for fixed px) is almost flat in the interval \y\ < ?7 cen , 
where \rj\ < i] ccn denotes the central rapidity region of the QGP. The expo- 
nent v is numerically found to increase with the initial temperature, hence 
increases with the total multiplicity rapidity distribution dN n /dy, which is an 
experimental observable. 

Thus, as the main conclusion, we propose that direct photons could potentially 
provide distinct experimental signatures of the transient nonequilibrium QGP 
created at RHIC and LHC energies, both in the form of a large enhancement 
at 1 — 2 < pt < 5 GeV/c as well as a power law transverse momentum 
distribution p^ v with an exponent v that is within the range 2.5 < v < 3 and 
increases with total multiplicity rapidity distribution dN n /dy. 

Discussions. 

(i) Our study of direct photon production focuses solely on the QGP phase and 
neglects contributions from pre-equilibrium stage as well as the mixed and 
hadronic phases, because we focus on a comparison between the usual equi- 
librium approach and the kinetic real-time approach. While the usual equi- 
librium approach treats the thermalized state as stationary, the real time 
kinetic approach allows to follow the time evolution of the density matrix 
as an initial value problem consistent with hydrodynamics. Furthermore, 
as discussed elsewhere in the literature, since most of the high-pj- photons 
originate from the very early, hot stage in the QGP phase [8,9,11,14] and 
photons produced in the mixed phase as well as the subsequent hadronic 
phase are mainly in the lower-pp region [18], the nonequilibrium yield from 
the QGP phase contributes dominantly to the total high-p T photons. 

(ii) While we have provided a systematic real-time description compatible with 
the initial value problem associated with hydrodynamic evolution, more 
needs to be understood for a complete description of all the different stages. 
The parton cascade approach to describe the early pre-equilibrium stage 
after the collision is an important first step in a full microscopic descrip- 
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tion, but perhaps a more consistent description of the initial condition for 
QGP formation must be based on the recent notions of a color glass conden- 
sate [35,36]. Hence, a complete treatment of the direct photon yield must, in 
principle, begin from the initial stage, possibly a color glass condensate, and 
obtain the real-time evolution of photon production. Clearly there must be 
many more advances in this field before such a program becomes feasible. 
The finite-temperature equilibrium calculations assume that the equilibrium 
thermal state always prevailed, thus ignores completely not only the initial 
stages but also the time evolution. Our approach while incorporating the 
time evolution during the stage of local thermodynamic equilibrium consis- 
tently, also neglects the initial stage. However, the advantage of our method 
is that if there is an estimate for the photon distribution at the onset of 
thermalization we could evolve the full quantum kinetic equation, including 
the higher order corrections in a as discussed in the text. We have focused 
on direct photons from a transient QGP for a direct comparison with equi- 
librium calculations, but obviously photons will continue to be produced 
during the mixed and hadronic phases. A detailed understanding of the 
phase transition as well as the hadronic photon production matrix elements 
is necessary for a more reliable estimate of the potential nonequilibrium 
effects after hadronization. 

(iii) Although we do not claim to have provided a completely detailed under- 
standing of potential nonequilibrium effects due to the lack of knowledge of 
initial conditions in heavy ion collisions, we do claim that our approach pro- 
vides a more systematic description of the finite-lifetime effects associated 
with a transient QGP. These effects lead to a distinct experimental predic- 
tion: a power law falloff of the distribution dN/d 2 pTdy near the central ra- 
pidity region which is distinguishable from a thermal tail for 1 — 2 < pt < 5 
GeV/c. We also note that the parton cascade results of Srivastava and 
Geiger [22] also seem to reveal a power law falloff of the photon distribution 
in this range of momentum with a similar amplitude and comparable expo- 
nent, as can be gleaned from Figs. 2 and 3 of Ref. [22]. Perhaps these two 
effects, namely the pre-equilibrium prompt photons and the direct photons 
from a hydrodynamically expanding QGP with a finite lifetime cannot be 
distinguished experimentally. Nevertheless, we emphasize that a power law 
departure from a thermal tail in the direct photon spectrum may very well 
be explained by nonequilibrium effects, either by a finite QGP lifetime as 
advocated in this article or by prompt photons from the pre-equilibrium 
stage. 

(iv) As many recent investigations [33,34] have suggested that the QGP pro- 
duced at RHIC and LHC energies is not expected to be in local chemical 
equilibrium, i.e., the distribution functions of quarks and gluons will prob- 
ably be undersaturated. An important extension of our work will be the 
study of nonequilibrium effects on direct photon production from a chemi- 
cally nonequilibrated QGP. 

(v) Transient QGP vs high energy particle collisions. An important and 
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very relevant question is that why not treat high energy particle collisions 
in the same manner, i.e., with the real-time evolution rather than with the 
S'-matrix approach. The answer to this question hinges on the issue of time 
scales. In a typical high energy particle collider experiment, the colliding 
"beams" are actually bunches or packets with a typical spatial extent of 
order 1 — 10 cm and hence the typical total time interval for the collision is 
of order 10~ 9 sec [37]. This time scale is many orders of magnitude longer 
than the typical hadronic interaction time scale ~ 1CT 23 sec, thus taking the 
infinite time limit is amply justified. This, of course, is the basis for using the 
S- matrix calculation in terms of asymptotic in and out states (at t = =Foo): 
the total interaction time (£/ — ti) is much, much longer than the typical 
hadronic time scale. With regard to the heavy ion collision, the consensus 
is that after an initial pre-equilibrium stage, a locally thermalized QGP 
is formed. It then evolves hydrodynamically, hadronizes and eventually the 
freeze-out of the hadronic gas ensues. Current theoretical estimates at RHIC 
energies indicate a total time between formation and freeze-out of order 
50 — 100 fm/c, with a QGP phase lasting for about 10 fm/c. These time 
scales, at least the lifetime of the QGP, are not several orders of magnitude 
larger than the hadronic interaction time scale, thus the infinite time limit 
taken in the S-matrix calculation is at best questionable. Therefore while in 
typical particle collider experiments the S'-matrix approach is valid (as has 
been demonstrated by over half a century of experiments!), the transient 
QGP with a finite lifetime of the order of the hadronic time scales merits a 
different treatment, which is the point of this article, 
(vi) An important and related aspect that requires further investigation is the 
finite size of the QGP. Much in the same way as the finite lifetime introduces 
nonequilibrium effects associated with energy nonconserving transitions, we 
expect that the finite size ~ 7 fm of the QGP will introduce uncertainty 
in the momentum of the emitted particles. This is clearly an important 
topic that deserves further and deeper study but is beyond the scope of this 
article. However, at this stage we speculate that if such effects are present 
they could bear an imprint in transverse flow. 

Work along many of these directions is currently in progress. 
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A Photon production rate and current correlation function 



The equilibrium thermal photon production rate is related to the Fourier trans- 
form of the thermal expectation value of the quark electromagnetic current 
correlation function [6,8,14]. In this appendix we show that the photon pro- 
duction rate given by Eq. (4.7) is related to the Fourier transform of the 
nonequilibrium expectation value of the quark electromagnetic current corre- 
lation function. 

In terms of the (Abelian) gauge invariant fields the interaction Lagrangian 
between quarks and photons reads [26] 

C iQt = eJ-A T , (A.l) 

— * — * 

with At the transverse component of the photon field and J the spatial com- 
ponent of the quark electromagnetic current 

J=J2e f ^ f ^ f . (A.2) 
/=i 

The invariant photon production rate given by Eq. (4.7) can be written in 
terms of J as 



E w§- X = itk? {i ~ iE ) ^ t] ■ ^ + c - c - (A - 3) 

It is convenient to introduce the Fourier transform of the transverse component 
of the quark current correlation function W^(u,p) defined by 



\ V ij (p) (j i+ (p, t)Ji~(-p, t')) = £j ^ W<(u,p) e-*> (t - t ' ) , (A.4) 

\ V ij (p) (.T-(p, t)Ji~(-p, t')) = £j ^ [W>(u,p)6(t' - t) 

+ W<(uj,p)6(t-t')] e --(*-*'), (A.5) 

where V^{p) = 5^ — p l p> is the transverse projector. A diagrammatic expan- 
sion [see Fig. 1(a)] shows that to order a and to all orders in a s Eq. (A. 3) can 
be expressed in terms of W < (u,p) as 



30 



where use has been made of the vacuum photon propagators. For a QGP with 
an initial state at time ij described by a thermal density matrix, one can show 
that the perturbative Kubo-Martin-Schwinger (KMS) condition holds to all 
orders in a s , i.e., 

W£(w,p) = e~" /% W£(u,p), (A.7) 

where Tj is the initial temperature. Consequently, one finds to all orders in 
a s that Wt(u,p) is a real quantity and related to the imaginary part of the 
retarded photon self-energy. Therefore Eq. (A. 6) becomes 

which is correct to order a and to all orders in a s . This result generalizes that 
of Refs. [6,8,14] to the nonequilibrium situation in the real-time formulation. 
Upon comparing Eq. (A.8) with Eq. (4.8), one finds 

K(u) = jW£(u,p) | onc . loop . (A.9) 
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